Light-induced valley currents and magnetization in graphene rings 
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We study the non- equilibrium dynamics in a mesoscopic graphene ring excited by picoseconds 
shaped electromagnetic pulses. We predict an ultrafast buildup of charge polarization, currents and 
orbital magnetization. Applying the light pulses identified here, non-equilibrium valley currents are 
generated in a graphene ring threaded by a stationary magnetic flux. We predict a finite graphene 
ring magnetization even for a vanishing charge current; the magnetization emerges due to the light- 
induced difference of the valley populations. 
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Introduction. - Since the recent fabrication of graphene, 
a monolayer of carbon, a number of fascinating phenom- 
ena have been uncovered, mostly owing to the quasi- 
rclativistic behavior of the carriers and their high mo- 
bility P, 0, S 0, S| ■ Two er-bonded interpenetrating tri- 
angular sublattices, A and B, build the graphene honey- 
comb lattice. The it and ir* bands govern the electronic 
properties near the neutrality point and result in coni- 
cal valleys touching at the high symmetry points K and 
K' of the Brillouin zone (BZ). Near K and K' the en- 
ergy dispersion is linear and the electronic properties are 
well described by the effective Dirac-Weyl Hamiltonian 
Hq @, H, IU [13] • The stationary states are degenerate in 
the spin S and the valley quantum numbers r = ±1. The 
latter correspond to the two non-equivalent if-points in 
BZ. Due to the suppressed intervalley scattering the con- 
trol of r eigenstates may be utilized for novel electronic 
0, 0] and optoelectronic applications Q. New physi- 
cal effects emerge due to confinement. E.g., in a meso- 
scopic graphene rings pierced by a magnetic flux, the ring 
confinement breaks the valley degeneracy and results in 
the persistent current Q . Experimentally, such graphene 
rings were fabricated and the Aharonov-Bohm effect was 
observed [loj ] . 

While a large body of work has been devoted to var- 
ious equilibrium electronic and optical properties, the 
non-equilibrium time-dependent phenomena in graphene 



are much less explored [111 1X21 ] . The present paper 



presents the first study on the non-equilibrium dynamics 
in graphene rings driven by asymmetric monocycle elec- 
tromagnetic pulses. Charge polarization and current car- 
rying states build up within picoseconds and are tunable 
by the parameters of the driving field. The states may 
become valley polarized resulting in a non-equilibrium 
valley currents. The valley population together with the 
charge current determine the magnetization of the ring. 

Stationary states.- We consider a graphene ring 0] of 
radius tq and width W (cf. Fig a)) threaded by a mag- 
netic flux of a strength <I>. As in [3,[lJ,[3|, the Dirac elec- 
trons are confined to the ring by the potential rV(r)a z 
at the boundaries, as resulting e.g. from a substrate po- 



tential [15|, [l6| • The polar-coordinates ring hamiltonian 
is [29] 
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+ rV(r)a z , (1) 



where o~ r = a • e r and er^ = a • e$ with e r and being 
the basis vectors of the polar coordinate system and o~ z 
is the Pauli matrix expressed in the pseudospin states of 
the two sublattices. v = 10 6 m/s is the Fermi velocity 
and $ = &/&q, where $o is the flux quantum. The 
eigenstates of H and J z (the z component of the total 
angular momentum with eigenvalues m) are 



i) m (r,4>) = R + {r)e^ m - 1 l^x + + R-{rV (m+1/2 ^X- , 

(2) 

where m = ±i,±|,... and a z X± — ±\x±- The gen- 
eral form for the radial parts is R+{r) = c\H^L > _i{r) + 

c 2 H^ ,(r) and R_(r) 

m 2 

where 
H 



(1)[(2)] 



= r\E\/v is a normalized radial coordinate, 
is the Hankel function of the first [second] kind, 
s = sgn(.E) selects the solution of the positive or the 
negative energy branch, and m = m + The bound- 
ary conditions and the normalization fix the coefficients 
ci and C2- For V(r) = if r 6 (r — ¥f,r + ^-), and 
V(r) — > +oo outside the ring [{| we find = ±T(T^ip at 
r = tq ± ^- . With this Eq. ^ can be solved numerically 
or for Wj ?'o <C 1 analytically [j| yielding the spectrum 



El 



SSyi ~\~ S\ 

Tiv 



m - 



cf> ST ™ 



sXr, 



1/2)2 
1 



,(3) 
■(4) 



n = 0,1,2,. 



and $ s c fi n = $ 



For fixed s, r, 



2 (n+l/2)ir ' 

and n the quantity modifies the energy spectrum 

as an effective (normalized) magnetic flux. The shift of 
the effective magnetic flux from $ has a different sign 
depending on the valley (r — ±1). For W/ro « 1 we 
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where 



r = (r — ro + tt) /W S (0, 1). For applications 
involving tunneling from the ring it is important to in- 
spect the case of a finite barrier boundary, i.e. V = Vq 
for r outside of the ring. To a first order of 7 = ^? <C 1 
we find that Eq. ((4]) applies with e n being replaced by 
(1 — 7)e„ and A„ by (1 + j)X n . For W = 0.1 /^m the 
condition 7 <C 1 means V ^S> 7 meV. Hence, our theory 
developed below is valid also for a finite barrier graphene 
ring. In a particular example of the boron nitride sub- 
strate we have V = 53 meV [ij} and therefore 7 = 0.13. 
Having specified the stationary single-particle states we 
proceed with the non-equilibrium calculations [3lj ]. 

Pulse-induced polarization.- To drive the non- 
equilibrium states in graphene rings we utilize asymmet- 
ric monocycle pulses, so-called half-cycle pulses (HCPs) 
3 III- For pulse duration t p shorter than the carriers 
characteristic time scale [32| the impulsive approxima- 
tion (IA) applies, meaning that the time-dependent car- 
rier wavefunction ^(r,(p;t) propagates stroboscopically 
as US 




FIG. 1: Dependence of the dipole moment generated in the 
graphene ring on the time past after the excitation and the 
normalized magnetic flux $ = 4>/'3?o G [ — 1/2, 1/2] (outside of 
this range the periodicity by $ — > $ + 1 can be used) in the 
case of 8 carriers in the ring at T = for (a) a — 1 and (b) 
a — 5. 
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where t~ and t + refer to times before and after the ap- 
plication of the pulse and a is the action delivered to a 
ring carrier of charge e by the HCP electric field F(t). 
The pulse triggers a time-dependent carrier density dis- 
tribution which depends on r = ±1, i.e. it is different 
for the two valleys. As a physical consequence, a time- 
dependent charge dipole moment is created in the ring. 
For the post-pulse dipole moment //™ (t) associated with 
a carrier startin g fr om the stationary state t = tq and 
m = toq we find [33| 



Mmn(^) = er o a h(Q) sin — cos 
to 
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where to = = Jo(^) + ^(- T ) de- 

notes the Bessel function of the order I, and £1 — 
ay/2 [1 — cos(2t/fo)J. The total electric dipole created 
in the ring for a fixed r, spin value iS, and t > is 
M r (*) = Em/>mW. where is the equilibrium dis- 
tribution function. For N T carriers in a given valley t at 
zero temperature T = we carried out the summation 
over to analytically. For an arbitrary even or odd N T we 



find respectively 
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Both expressions apply for l>^ ff = $ - J e [-1/2, 1/2], 
outside of this interval n ST (t) is determined from the peri- 
odicity in $^ ff with a period 1 . The total dipole moment 
fi(t) depends on the distribution of the carriers between 
the valleys that in turn depends on the magnetic flux 
The spin degeneracy is also important. One can show 
that jumps in the population of particular states take 
place only at the points $ = — | j — — j —T, 0, r, 
in the interval [—1/2, 1/2], where we denote r 
The dynamics of the dipole moment for N = 8 carriers 
is shown in Fig. Q] as a function of the applied station- 
ary magnetic flux for two different excitation strengths 
a = 1 and a = 5, showing that the ring electric dipole 
and hence the associated light emission are controllable 
by <E> and a. For ro = 1 /im and HCPs with a sine-square 
shape and a time duration of 0.5 ps, a — 1 corresponds 
to the peak value of electric field F = 26 V/cm. 
Note, the boundary conditions break the effective time- 
reversal symmetry [9( making the states corresponding to 
different r but otherwise to the same quantum numbers 
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non-degenerate. The dynamics of the charge polarization 
for confined carriers is however the same in both r valleys 
for $ = 0. This follows from the invariance of the states 
under r — > — r, and m — > —to at $ = 0, as evidenced 
by Eqs. ©. This degeneracy is lifted by applying a sta- 
tionary magnetic flux $ ^ 0. The density distribution of 
carriers in the ring becomes valley-polarized. 
N on- equilibrium charge and valley currents.- The electric 
current density j = evip^aip has the ^-component of the 
current density 



U(r) = -21m[R + (r)R*_(r)] 



(10) 



The total charge current is I = J j^dr. To a zero or- 
der in W/ro « 1 we find I — 0. Only higher order 
corrections in W/ro give rise to a non- vanishing ring 
current. For the eigenstate specified by s,T,n,m the 
lowest order correction 0] to the current follows from 
Inm — ~dEntp I ^ using the energy spectrum in the con- 
sidered limit [34|. For n = 0, s = 1 this current is equal 



to II = -I 
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L - L - 5 — . I 1 or a ring 



with ro = 1 /im and W — 100 nm we obtain Iq = 0.16 nA, 
if ro = 425 nm and W = 150 nm as in Ref. llCj we find 
Iq = 1 nA. In both cases I A is valid if t p <C 3 ps. Such 
HCPs are experimentally available The total cur- 

rent in the ring I is the sum of an equilibrium (persistent) 
current J cq and a non-equilibrium time-dependent cur- 
rent part I neq (t) generated in the ring: I = J cq + I neq (t). 
The equilibrium part is given by I eq = J2mr fmr^m- F° r 
T = it is given in Ref. H for N = 1, 2, 3, 4. We derive it 
for any TV in n = 0. 

A non-equilibrium ring current is generated by a 
sequence of two time-delayed mutually perpendicular 
HCPs (see Fig.^a)), similarly to the pulse-current gener- 
ation in semiconductor rings SESii- This scheme 
allows for shorter excitation times compared to the res- 
onant excitation schemes using circular polarized pulses 
At t 



27] 



we apply linearly polarized 
(along the x-axis) pulse that creates a time-dependent 
charge polarization along the x-axis (cf. Fig. [IJ. The 
second pulse is linearly polarized along the y-axis and is 
applied at t = t y . It generates a non-equilibrium current 
depending on the charge polarization created by the first 
HCP. The delay time should be short enough so that re- 
laxation processes are negligible in between the pulses. 
In the IA The generated non-equilibrium current reads 



J„ 



er 



-i e(t-t y ), 



(11) 



where a y is the excitation strength of the second HCP 
and H x {ty) is the dipole moment created by the first HCP 
just before the application of the second HCP. Equation 
(fTTjl delivers the total current as well as the individual 
currents in each of the two valleys, in which case fJ>x(ty) 
should be associated with the charge carriers in the re- 
spective valley. Defining the valley current as the differ- 
ence between the currents flowing in two opposite valleys 



divided by the particle charge we find for the generated 
valley current 



ncq 



er 



e(t-t v ). 



(12) 



On a longer time scale set by the relaxation processes 
the non-equilibrium current decays due to dissipation. 
Thereby the incoherent electron- pho non scattering plays 
usually the most important role [23j, |24|. Specifically for 
a free-standing graphene, scattering by flexural phonons 
is dominant at low temperatures [281 ] . 
An example of the dependence of the generated total 
charge current on the delay time t y is depicted in the up- 
per panel of Fig. [3^b). The oscillating character of this 
dependence is determined by the dynamics of the dipole 
moment generated by the first HCP. The lower panel of 
Fig. [Hb) demonstrates the dependence of the generated 
valley current on the delay time t y . This current arises 
as a consequence of the different contributions to the to- 
tal dipole moment from the two different valleys in pres- 
ence of a static magnetic flux (here we used $ = I/tt) 
at the time moment t — t y . Comparing the upper and 
the lower panels we conclude that tuning the pulses de- 
lay may result in a vanishing total generated current I neq 
while the generated valley current 7„ eq is finite. It is also 
possible to create I n cq 7^ with I^ cq = 0. Under the 
conditions of Fig. EJJb) the generated currents have the 
same order of magnitude as the persistent currents. The 
non-equilibrium contributions are enhanced however by 
increasing the HCPs excitation strengths. An increase of 
the excitation strength a x of the first HCP beyond the 
values around 1 does not lead however to an increase of 



I (T 



H under the conditions where 7^ eq 



(In 



ishes because for this, certain delay times are required. 
In the strong excitation regime the nonlinear oscillations 
of the dipole moment collapse [23| shortly after the exci- 
tation (cf. Figs. QJa) and (b) in the range t/to € [0,2]). 
For a further increase of the currents under these condi- 
tions a y should be increased. 

The ring charge current is associated with a mag- 
netic dipole moment via M = \ J f x j d 2 f, i.e. M = 
ir J jj,r 2 dr. From Eqs. (jTOJ) , (J5j) and ((6|) we infer for the 
non- vanishing lowest order of W/ro 
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Tir Q I 
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(13) 



where Q sn — — is the difference in the valley 
population for fixed s and n. For a vanishing total cur- 
rent in the ring and s = 1, n = 0, Eq. (|13| simplifies to 
M = 4:Qirr 2 I . Note, the valley polarized magnetic mo- 
ment is also a generic feature of the monolayer graphene 
with a broken inversion symmetry (e.g. due to the action 
of the substrate potential) 0]. The difference in the val- 
ley population in Eq. (TT3|) arises in equilibrium for certain 
ranges of <I> ^ 0. It can be also generated e.g. by injec- 
tion of external non-equilibrium carriers to the graphene 
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(a) 5 (b) 



FIG. 2: (a) Current generation in the ring by application of 
two HCPs polarized along mutually perpendicular directions 
and delayed in respect to each other by the time t y . (b) Up- 
per and lower figures show the total current and the valley 
current, respectively, generated in the graphene ring in de- 
pendence on the delay time t y . Value of the magnetic flux 
is set to $ = $o/t, number of carriers is N = 8, excita- 
tion strengths of both HCPs are equal to a = 1. Arrows at 
t y = to indicates the delay time for which a valley current 
with no total charge current is generated whereas for a delay 
time t y = 2f^io a total charge current with equal contribu- 
tions from both valleys is generated. 

ring, opening thus a new way for an ultrafast detection 
of the valley number. Finally, we note our results are 
valid for weak pulses in which case a small angular popu- 
lation around the ground state is created and many-body 
effects remain subsidiary. Strong excitations go beyond 
the present model and the influence of many-body inter- 
actions may decisively alter the above predictions. 
Conclusion.- Short linearly polarized asymmetric light 
pulses trigger a non-equilibrium carrier dynamics in 
graphene rings threaded by a magnetic flux. The induced 
charge polarization is detectable by monitoring the emit- 
ted radiation. Delayed pulses with different polarization 
axes drive non-equilibrium charge currents and hence an 
orbital magnetization. For appropriate pulses, equal con- 
tributions from both valleys is achievable as well as pure 
valley currents. The ring magnetization depends on the 
difference in the valley population. The predicted effect 
is operational in presence of tunneling allowing thus for 
swift injection or detection (via ring magnetization) of 
valley currents in coupled graphene structures, e.g. wires, 
offering new realization of ultrafast valleytronics devices. 
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